The role of body rotation in bacterial flagellar bundling 
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In bacterial chemotaxis, E. coli cells drift up chemical gradients by a series of runs and tumbles. 
Runs are periods of directed swimming, and tumbles are abrupt changes in swimming direction. 
Near the beginning of each run, the rotating helical flagellar filaments which propel the cell form a 
bundle. Using resistive-force theory, we show that the counter-rotation of the cell body necessary for 
torque balance is sufficient to wrap the filaments into a bundle, even in the absence of the swirling 
flows produced by each individual filament. 

PACS numbers: 87.16.Qp, 87.16.-b, 46.70.Hg 



Although bacteria are among the simplest systems for 
the study of cell motility, many puzzles remain. Chief 
among these is the mechanics of the bundling and un- 
bundling of flagellar filaments in the chemotaxis behav- 
ior of bacteria such as E. coli and Salmonella ^, ||. 
These cells move toward higher concentrations of favor- 
able chemicals by executing a series of runs and tum- 
bles H . The runs are periods of directed swimming. At 
the end of each run, the cell randomizes its direction by 
tumbling. If the cell happens to head in a favorable direc- 
tion, the likelihood of tumbling reduces, making runs in 
this direction longer on average compared to runs in the 
unfavorable direction. Propulsion during a run is gener- 
ated by the rotation of several helical propellers, known 
as flagellar filaments. Unlike eukaryotic flagella bac- 
terial flagellar filaments are passive elements driven by 
rotary motors embedded in the cell wall. Near the be- 
ginning of a run, the motors turn in a counter-clockwise 
direction (when viewed from the outside of the cell) , and 
the left-handed filaments come together to form a bundle. 
At the end of a run, one or more of the motors reverses, 
and one or more of the filaments fly out of the bundle and 
cause the cell to tumble. This process is complex and in- 
volves changes in filament handedness and pitch. The cell 
soon sets out on a new course but regains its initial speed 
only after the aberrant motors have reversed again and 
their filaments have regained their normal conformation 
and rejoined the bundle W. 

Although qualitative partial explanations for bundle 
formation have appeared in the literature [| a 
mathematical theory has not. In this note we begin to 
construct this theory with a quantitative treatment of 
one aspect of the bundling phenomenon: the role of the 
cell body rotation. This rotation and the accompanying 
hydrodynamic resistance arise to balance the torque ex- 
erted by the rotating bundle on the cell body. Thus, in 
the body- fixed frame, there are two kinds of flows which 
contribute to bundling: the flow due to frame rotation, 
and the swirling flows set up by each individual filament. 
Here we focus on the flow due to rotation of the body- 



*Electronic address: Thomas Powcrs@brown.edu 



FIG. 1: Model problem. A naturally straight but flexible 
rod, initially parallel to the z-axis, with one end held fixed 
with zero moment at x — b, y — 0, z — 0. We seek the 
steady-state shape when the end of the rod is forced to rotate 
in the x-y plane about the z-axis at frequency ub . 



fixed frame; the swirling flows and interactions among 
flagellar filaments will be treated in a separate publica- 
tion H . Our treatment is in the spirit of Machin || , who 
used a similar approach to argue that eukaryotic flagella 
could not be passive elements driven by motors at the 
cell body (see also ref. |Hj). 

Fig. [l| illustrates the model problem. For simplicity, re- 
place the helices with straight but flexible rods of length 
L, rotated with the frequency lub about the cell body 
axis of symmetry, z. Let b denote the distance between 
the axis of the unstressed rod and the z-axis. Since the 
body is about a micron across, and flagellar filaments are 
typically six to ten microns long, we suppose b <C L. We 
also disregard the rotational disturbance flow arising from 
the no-slip condition at the cell body. In the body-fixed 
frame, this disturbance flow reduces the net rotational 
flow near the body. We disregard this disturbance flow 
since the flow field of a sphere of radius a rotating with 
angular velocity uj takes the form v = u>xr(a/r) 3 , falling 
off rapidly with r Likewise, it is argued below that 
the axial drag on a filament due to the nonzero swim- 
ming velocity plays little role in our problem. Finally, 
we focus our attention on the contribution to flagellar fil- 
ament wrapping due to body rotation, and not the flows 
set up by the individual rotating filaments, by ignoring 
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the hydrodynamic interactions among the rods. Thus, it 
suffices to consider the shape of a single rod. 

During runs, the left-handed flagellar filaments turn 
counter-clockwise (when viewed from outside the cell), 
and the body turns clockwise (when viewed from 
behind — i.e. from the distal end of the bundle). When 
our model filament is turned about the body-rotation 
axis z in this same sense (clockwise when viewed from 
the positive z-axis, see fig. |l|), it forms a right-handed 
shape (e.g. see fig. || and note that the proximal end 
x/L = b = 0.1 is in the plane z — 0, and the distal 
end with x near has positive z coordinate). Further- 
more, two lefthanded helices rotating about their respec- 
tive axes with proximal ends held stationary will lead 
to a flow which also tends to wrap the helices around 
each other in a righthanded manner @. Thus, the body- 
rotation effect treated here and the swirling flow effect 
treated in [|| act to wrap the filaments in the same sense. 

Since typical Reynolds numbers for swimming bacte- 
ria are of order 10~ 6 inertia is unimportant and the 
steady-state rod shape is determined by a balance of vis- 
cous and elastic forces per unit length. For gentle distor- 
tions of a slender body, the viscous forces per unit length 
are well-approximated by the resistive-force coefficients 
(e.g see and references therein): / = (±u±_ + C\\ u \\i 
where itj_ and U\\ are the perpendicular and parallel com- 
ponents of the local rod velocity relative to the fluid ve- 
locity v. u = dr/dt — v. The transverse friction coeffi- 
cient (per unit length) is of the form 
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log(L/o) + l/2' 



(1) 



where r\ is the fluid viscosity, a is the rod radius, and log 
denotes the natural logarithm . As discussed below, 
does not enter the analysis since we work in the linearized 
approximation. Resistive-force theory gives an accurate 
value for the drag force per unit length on a slender fil- 
ament except near the filament ends; however, the effect 
of this error on the shape is 0(a/L) @. There is also 
a viscous torque tending to twist the rod; however, the 
effects of this torque are subleading compared to the ef- 
fects of the translational drag |l3| . To see why, note that 
the total torque from rotational drag is 0{ujC, r L), where 
Q r = Aty^o 2 is the friction coefficient for rotation p[ . 
The total torque from translational drag is O(b(ojb(±)£), 
where, as we shall see below, only the portion of the rod 
within a distance I of the held end (z = 0) contributes 
to the translational drag. The ratio of these two torques 
is 0((a/b) 2 (L/e)[\og(L/a) + l/2}). For the representative 
values L = 10 /xm, a w 10 nm, and b ss 1 /^m, this ratio is 
small even if L/£ « 10. We therefore disregard rotational 
drag and twist strain. 

To find the bending force per unit length, note that 
since b <C L, the displacement of any rod element will also 
be small. Thus, the elastic energy is well-approximated 
by the quadratic expression 
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where x and y are as in fig. |], and A is the bend mod- 
ulus [fl5|| . Since the variation in rod shape is rapid for 
sufficiently high rotation rate, even for small b, this ap- 
proximation eventually fails and must be replaced by the 
full geometrically nonlinear elastic rod energy. As we dis- 
cuss below, the rotation rates of interest are small enough 
for (||) to hold. The variational derivative of (||) yields 
the elastic bending force per unit length: —8£/5r± = 
—Ad 4 r±/dz 4 , where r±=xx + yy. 

To leading order for b <C L, the motion of the rod is 
purely perpendicular to the rod centerline, yielding the 
equation of motion [pi 
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where v± is the transverse fluid velocity. Since inertia 
is unimportant in the limit of zero Reynolds number, 
eqn. (Q) applies equally well to the rotating frame in 
which the rod is fixed and the flow is v± = lo-qz x r±. 
Such a flow tends to wrap the rod around the z-axis in 
a shape with a helical modulation and exponential enve- 
lope. In the steady-state, eqn. (||) reduces to 
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where I = (A/^lub) 1 ^ 4 is the characteristic length scale 
associated with bending and drag ||. The solution to 
eqns. (|4|j^) is a simple generalization of Machin's solution 
to the in-plane bending problem: 



x(z) = ^A n exp(r„z/f), 



(6) 



where ri,...,rg are the eight eighth roots of —1. The 
wavelengths A„ and decay lengths v n of the eight fun- 
damental complex solutions exp(r n z/£) (with r n — 
2iri/X n + l/v n ) are comprised of the four possible com- 
binations of A„ = ±16.419 and v n = ±1.0824, and 
the four possible combinations of X n = ±6.8009 and 
v n = ±2.6131. 

The boundary conditions determine the amplitudes 
and phases of the coefficients A n . At the distal end 
z = L, there is zero force and moment: Ad 3 r±/dz 3 = 0, 
Ad 2 r±/dz 2 = At the proximal end, flagellar fila- 
ments are connected to the rotary motor by a hook which 
is more flexible than the rest of the filament. We simply 
model this flexible connection as a hinge with zero mo- 
ment at z = 0: Ad 2 r±/dz 2 = 0. (The other extreme, a 
rigid hook with dr±/dz = leads to qualitatively sim- 
ilar shapes for £/L < 1, except near z = 0.) Finally, 
r±(z = 0) = bx. Applying these boundary conditions 
to the solutions in eqn. (g) with b = L/10 yields the 
shapes shown in figs. |^-^. For large £, the rod is very 
stiff and does not bend; it is easy to show that in the 
limit oft/L > 1 that x(z) = 6(1 -3z/2) + 0((L/£) 4 ) and 
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FIG. 2: Projection of shapes of a rotating flexible rod onto 
the x-y plane for i/L = 0.1 (solid line), i/L = 0.2 (dashed 
line), and i/L = 0.5 (dot-dashed line). 



y(z) = 0((L/£) 4 ) for the hinged (zero moment) bound- 
ary condition at z = 0. In the lab frame, the rod pivots 
about the point z/L = 2/3, tracing out a cone (a rigid 
rod confined to the plane also pivots about this point in a 
viscous fluid pi). When i/L < 1, the rod spirals around 
the axis of rotation, with the spiral becoming more com- 
plete as £ gets smaller and smaller. Note the anisotropy; 
the projection of the shape x-y plane is elongated along 
the x-axis (fig. ||). The rod configuration is a compromise 
between minimizing the bending energy and minimizing 
the dissipation rate. For i/L 3> 1, elasticity dominates, 
and the rod is straight. For i/L <C 1, viscous effects 
dominate and the rod bends to align along the axis of 
rotation and thus minimize the dissipation rate. In this 
limit, the linear approximation for the shape of the rod 
becomes invalid, and we must replace (|^) with the full 
expression for the curvature. This nonlinear problem is 
readily solved by standard methods (see e.g. p6t); the re- 
sult is that inaccuracies of only a few percent arise when 

To assess the importance of the role of body rotation in 
bundling, we estimate the characteristic length i. Var- 
ious estimates have appeared for the flagellar filament 
stiffness A, from 1CT 24 N m 2 @ to 10" 22 N m 2 §. 
Fortunately, the characteristic length I is not very sensi- 
tive to the value of A. To estimate the transverse drag 
coefficient cc l n - we use the viscosity of water 
i] = 0.001 N s/m 2 , a typical length L = 10 /im, and 
a diameter 2a = 20 nm. With a typical body rotation 
rate of wb = 10 Hz fl and the range of stiffnesses quoted 
above, the characteristic length i is found to be two to six 
microns. Therefore, the filaments are sufficiently flexible 
for the observed body rotation rate to contribute signifi- 
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FIG. 3: Shapes of a rotating flexible rod, projected onto the 
z-x plane, for i/L = 0.1 (solid line), i/L — 0.2 (dashed line), 
and i/L = 0.5 (dot-dashed line). Vertical amplitudes have 
been exaggerated for clarity. 




FIG. 4: Shapes of a rotating flexible rod, projected onto the 
y-z plane, for i/L = 0.1 (solid line), i/L — 0.2 (dashed line), 
and i/L = 0.5 (dot-dashed line). Vertical amplitudes have 
been exaggerated for clarity. 

cantly to bundling. Furthermore, the linear treatment of 
the rod shape is justified. Presumably, body rotation is 
especially important for the bundles which include many 
right-handed filaments and a single left-handed filament, 
as observed in ref . Q . 

Including axial drag does not alter the conclusions. 
Axial drag due to the swimming velocity leads to a ten- 
sion gradient in the rod which slightly increases the spiral 
pitch. Assuming a constant tension equal to the maxi- 
mum tension at the base of the rod and disregarding the 
shadow effect of the cell body yields an upper bound on 
the change in pitch. For a swimming velocity of about 
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30 /^m/sec, the change in pitch is small compared to the 
pitch. 

The purpose of this work has been to point out the im- 
portance of body rotation for flagellar filament bundling. 
In order to focus on the essential physics of this element 
of the bundling phenomenon, we have disregarded several 
important but complementary effects, such as the helical 
shape of the flagellar filament and the flows induced by 
the individual filaments ||. Despite these simplifications, 
we have shown that bacterial flagellar filaments are flex- 



ible enough for body rotation to lead to wrapping. 
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